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(N ■ ABSTRACT 

r~| ■ It is shown that electric current flowing along an axially non-uniform magnetic flux tube produces an 
Oh associated non-linear, non-conservative axial MHD force which pumps plasma from regions where the flux 
' tube diameter is small to regions where it is large. In particular, this force will ingest plasma into the ends 
c/2 ■ of a fat, initially potential flux tube and then pump the ingested plasma towards the middle bulge, thereby 
^ causing mass accumulation at the bulge. The ingested plasma convects frozen-in toroidal magnetic flux 
Qh which accumulates at the middle as well. Flux accumulation at the bulge has the remarkable consequence 
CO ■ of causing the bulge to diminish so that the flux tube becomes axially uniform as observed in coronal loops. 
. i-H Stagnation of the convergent plasma flow at the middle heats the plasma. A small number of tail particles 
bouncing synchronously between approaching fluid elements can be Fermi-accelerated to very high energies. 
' Since driving a current along a flux tube is tantamount to helicity injection into the flux tube, this mass 
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ingestion, heating, and straightening should be ubiquitous to helicity injection processes. 
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fSJ . Remarkably detailed images of coronal loops provided by the TRACE spacecraft [Aschwanden et al 
(2000)] indicate that these loops have cross-sectional area varying by only 10-20% over their entire length; 
this surprising behavior cannot be explained by potential or force-free field models [Klimchuk (2000]. We 
present here a model explaining this behavior. Since the loops are typically twisted by a fraction of a turn 



CJ , over their length, they contain magnetic helicity and so our model should be intrinsic to helicity injection 
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cor processes 



A coronal loop will be represented here by an axisymmetric flux tube with straight axis (cf. top, Fig. 1) 
O and cylindrical coordinates (r, <p, z) will be used. The loop is characterized by an initially potential poloidal 
flux function 

ip{r,z)= / B z (r',z)2irr'dr' (1) 



& ' with associated poloidal field 

B P oi = x V(f>. (2) 

Axial non-uniformity corresponds to ip being z-dependent and bulging corresponds to ip~ 1 d 2 'ip/dz 2 > 0. 
Similarly, the poloidal current is written as 



I(r,z) = I J z {r',z)2Txr'dr' (3) 
with associated poloidal current density 



o 



Jpol = TT^I X V<t>. (4) 
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We show that establishing a steady-state current / involves three sequential stages having distinct 
physics. The first stage, "ramp-up", has physics akin to a linear Alfven wave, but it is assumed that 
the ramp- up rate is sufficiently slow that the effective Alfven wavelength is infinite. This means that 
retarded time effects due to wave propagation issues are negligible and the current ramps up everywhere 
simultaneously as in an ordinary electrical circuit. The second stage, "axial flow", has dl/dt = 0, but is 
not in MHD equilibrium because unbalanced, non-conservative J x B forces exist which drive plasma flows. 
The third stage, "stagnation", involves convection of magnetic flux by the flows, plasma heating as a result 
of flow stagnation, and straightening of the ip profile until MHD equilibrium is established. 



First stage (Ramp-up) 

We represent the current ramp-up by the time-dependence 

1 + tanh(t/r) 

I[t) = !o ^ (5) 

where the ramp-up time is assumed to be much longer than the time it takes for an Alfven wave to propagate 
the length h of the flux tube, i.e., r >> h/vA- Alfven wave propagation effects are therefore unimportant 
in which case the system behaves like an electric circuit. From Ampere's law the toroidal magnetic field is 

p I(r,z,t) 

B d r i z it) = — ^ — • (6) 

The toroidal component of Faraday's law is 

dE r dE z dB A 



dz dr dt 



(7) 



We note that B A has minimal z dependence and that Ohm's law implies E z ~ 0. Thus, integration of Eq.(0) 
with respect to z gives 

Er ^~ 2^ at (8) 

where on the basis of symmetry the location z = is set to be at the axial midpoint of the flux tube. Since 
there is no axial force in this stage, U z remains zero and so the radial component of the ideal Ohm's law 
gives 

U<P = ~B L Z ~ 1 h^B~ z W (9) 

showing that Us is finite only when I is changing. The change in Uj, (toroidal acceleration) implies the 
existence of a radial current determined from the toroidal component of the equation of motion 

= ~ JrBz - (10) 

This current is just the polarization current [Chen(1984), p. 40] 

j_ p dU^ _ p 8E r pzp d 2 I 

r ~ B z 8t B 2 Z 8t ~ 2nrBl dt 2 ' [ ' 

The transient toroidal velocity given by Eq. (||) results in an azimuthal displacement of the plasma, 

showing that the plasma motion in this stage follows the twisting of the magnetic field (the field line can 
be thought of as being frozen to the plasma so that when the field line twists, so does the plasma). Thus 
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Initially untwisted potential flux loop 




Fig. 1. Top: Initially potential flux tube. Bottom: Flux tube with axial current has same poloidal profile, but is 
twisted. Cylindrical geometry (r,<f>,z) is used; z = is the flux tube axial midpoint and the flux tube ends are 
at z = ±h. 



J r is finite only when I is changing and J r is first negative and then positive, corresponding to toroidal 
acceleration followed by toroidal deceleration. The r direction here is really a proxy for the Vip direction, 
since the polarization current is in the Vip direction (because the flux tube is long and slender, the r direction 
is nearly the same as the Vip direction). 

Once / has been established, both Ua, and the polarization current J r remain zero. The poloidal flux 
function if; at this stage has not changed from its original value — all that has happened is that a toroidal 
field has been added so that the total magnetic field is now 

B4WxV f i|/ 1 oIVi (13) 

The original potential flux tube has become twisted as shown in the bottom sketch of Fig. 1. Since if) is 
unchanged from its potential value, the bottom sketch in Fig.l has the same poloidal profile (envelope) as 
the top sketch. This means that 

(Vib x Vd>) r ( 1 \ 
= rV<f> ■ V x 1 \^ 91 = --V .(-Vip)=0. (14) 
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Second stage (flow acceleration) 

The magnetic force is 

" "iter J 



X B= (J po i + 3 tar) X (B po i + B 
= Jpoi x -Bp G ; + Jp l X Btor 



'tor + Jtor X B po ;. (15) 



The term J po ;xB po ; provides acceleration in the toroidal direction and involves the component of J po i which 
is perpendicular to B po z. However, in the previous section we showed that the component of J po i which is 
perpendicular to H po i is the polarization current and that this current scales as d 2 I/dt 2 . Thus, when / is in 
steady state there is no polarization current and no toroidal acceleration, and so U<f, remains zero. It was 
also shown that Jtor = J<f>$ is zero so long as the poloidal flux surfaces are unperturbed from their initial 
potential shape. We assume in this second stage that (i) the current is constant in which case J po ;xB po ; = 
and 17^ = and (ii) the poloidal flux surfaces are unperturbed from their initial potential shape in which 
case Jtor = 0. Thus, during this second stage the magnetic force reduces to 



J x B = J pol x B tor = -L (VV x V0) x = -t-^V/ 2 



(16) 



The above force is non-conservative (i.e., V x (J x B) is non-zero) and so cannot be balanced by a pressure 
gradient since a pressure gradient is conservative (i.e., V x VP = 0). Thus, it is not possible for equilibrium 
to be attained in this stage. The only way for an equilibrium to be achieved is for the poloidal profile of the 
magnetic field to change, which is what happens in the third stage, to be discussed later. 

The fact that J po ;xB po ; = means that (VI x V</>) x (Vip x V0) = which in turn implies that V/ is 
parallel to and so I must be a function of vp, i.e., I = I{ip). Thus, the poloidal current flows along the 
poloidal flux surfaces. This is consistent with the well-known Hamiltonian dynamics concept that, because of 
conservation of canonical angular momentum, particles in a toroidally symmetric geometry cannot make an 
excursion exceeding a poloidal Larmor radius from a poloidal flux surface [e.g., see p. 207-208 of Bellan(2000)]. 
In other words, Hamiltonian mechanics forbids the existence of steady current in the direction normal to a 
poloidal flux surface. 

As sketched in Fig. 1, the poloidal flux function is bulged near z = 0, corresponding to a weaker magnetic 
field near z = than at the ends z = ±h. This would be the situation if the source currents for the poloidal 
field were located external to the flux tube and so the middle of the flux tube would be further from the 
source currents than the ends. Since I = I (if)), the current channel would also be bulged. 

Equation (|l6|) implies that the z component of the equation of motion is 

dU z ,, „, 1 dl 2 d ( B 2 \ 

= ( J x B '- = (17> 

This means that there is a force accelerating plasma from regions where is strong to regions where i?| 
is weak. Since I = I(ip) and since ip is bulged in the middle, must be stronger near z = ±h where 
the current channel diameter is small than at z = where the current channel diameter is large. There 
consequently must be an acceleration of plasma from both ends (i.e., z = ±/t) towards the middle (i.e., z = 0) 
as shown in Fig. 2. The convergent axial pumping is similar to the "sweeping magnetic twist mechanism" 
discussed by Uchida and Shibata (1988), but it should not be considered a wave because it involves actual 
convection of material and not propagation of energy through a material. We note in passing that there 
could be a few exceptional particles collisionally bouncing back and forth between the approaching fluid 
elements. These exceptional particles would be accelerated to very high energy by the Fermi acceleration 
process, and so one would expect to see a tail of energetic particles develop in the vicinity of z = 0. The 
Fermi process would thus predict that the most energetic particles would be located around the top of an 
arched coronal loop and such is indeed what is observed (Feldman, 2002). 
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Third stage (stagnation, heating, and straightening) 

The flows from both ends converge at the middle and must come to a halt at z = 0. Convergence of 
flows means that V • U is negative. This has important implications for the magnetic field as can be seen 
by considering the induction equation toroidal component, 

^ = r B por V (?f) - rlVV - B^ ■ U pol . (18) 

We have shown that (i) U<j, = 0, (ii) U po ; — > at the stagnation layer, and (iii) V • XJ po i is negative. Thus, 
in the vicinity of the stagnation layer 

= -B^V ■ XJ pol (19) 

showing that must increase at the stagnation layer (increase of magnetic field at regions of local flow 
convergence has been discussed in a more general context by Polygiannakis and Moussas (1999)). 

The continuity equation in the vicinity of the stagnation layer gives V • XJ po i = —p~ 1 dp/dt and so in the 
vicinity of the stagnation layer 

1 dB 4> = l^P m) 
Bff, dt p dt' { ' 

Thus, in the vicinity of the stagnation layer increases in proportion to the accumulation of mass. Since I is 
constant during this stage, if Ba, increases, the radius of the current channel must decrease so as to maintain 
PqI = 2irrB^ constant. However, since I = I(tp), if the radius of the current channel decreases, then the 
radius of tp must also decrease, thereby reducing the bulging. Ultimately, the bulge becomes vanishingly 
small as more and more mass accumulates at the stagnation point and eventually a plasma-loaded, axially 
uniform flux tube results. 

Since stagnation involves conversion of flow velocity into thermal energy, the plasma will be heated at 
the stagnation layer. The resulting plasma temperature can be estimated from the details of the stagnation 
process. The axial magnetic force in Eq.([l7]) involves the gradient of the toroidal field energy density and 
so the kinetic energy gained by the plasma will be of the order of the change in toroidal field energy density 
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B 2 /2fiQ between the ends [z = ±h) and the midpoint z = 0. Because B 2 /2po is larger at the ends, the 
plasma flow kinetic energy is pU 2 /2 ~ P?/2/Uo and at the stagnation layer z = this flow energy is 

L ™ J z=±h 

converted into heat so pU 2 = B 2 /po — > nvn e v\ e + nrriiV^. 

Since the poloidal field is much stronger than the toroidal field and is approximately B z , this means 
that the flow stagnation causes the plasma to develop a state where 

2p (nm e v^ e + nm^-) _ -g| / ^q! \ 2 _ / p I\ 2 a 2 



where a is the radius of the current channel. However, pqI = J ds ■ V x B and rp = J ds ■ B where the 
surface integral is over the cross-sectional area of the flux tube. If we define the ratio of poloidal current to 
poloidal flux as 

a = fiol/t/j, (22) 

the (3 predicted from flow stagnation is 

0= aV/2. (23) 

Thus, when a current is made to flow along an initially bulging current-free current channel, the 
current channel will twist up (helicity injection), plasma will be ingested from both ends, accelerated toward 
the middle where it accumulates and heats up at a stagnation layer. The equilibrium will become straight 
(filamentary) and have (3 = a 2 a 2 /2 where a = hqI/i/j and a is the radius of the current channel. To an 
outside observer the current channel will look field-aligned since the current is axially uniform and appears 
to be embedded in an axially uniform axial magnetic field. However, oblivious to the outside observer, within 
the current channel there is a Bennett pinch-like radial force balance between plasma pressure pushing out 
and JzBfp magnetic force pushing in. 

Although the predicted (3 is typically very small, its effect is crucial. To see this, consider that 
in equilibrium J x B = VP so that B-VP = (2-it)^ 1 (Vip xVcp + HoIV(f>) -VP = and so P must be a 
function of ■0, i.e., P = P{i[>). Defining i/jq as the flux on the flux surface where P vanishes, we can write 
P(r,z) = (1 — t/j(r, z)/ipo)Po where Po is the on-axis pressure (i.e., where if) = 0). We can also write 
PoI( r i z ) = cttp(r, z) and so J x B = VP can be written in Grad-Shafranov form as 

d (ld^\ d 2 V> a 2 2 P ( 

r TT + ~o~T + a ^ = 47r w ~r- ( 24 ) 

or \r or J oz z ipo 

If a 2 = Att 2 fiQa 2 Pq / tpQ where a is the flux tube radius at z = 0, then the only solution to Eq.(^3|) satisfying the 
specified boundary condition that P vanishes when ip = i/jq is the particular solution ip(r, z) = tpor 2 /a 2 . This 
means that the flux tube must be axially uniform when a 2 = 47r 2 poa 2 Po / V>q . Defining Bq = ipo/ira 2 as the 
axial field at z = 0, it is seen that this condition for axial uniformity corresponds to a 2 = Att 2 p^a 2 Pq / B^ir 2 a A 
or a 2 a 2 /2 = j3 where j3 = 2/iqPq/ Bq. This equilibrium has = 0, so all confinement is provided by the 
Bennett pinch force ~ JzB^. The current is purely in the z direction, but the magnetic field is helical. 

The situation of small but finite (3 is substantially different from the case of zero (3 because the system 
is constrained to be axially uniform if and only if (3 = a 2 a 2 /2. The arguments presented in the discussion 



of Eqs.(21-23) show that the MHD dynamical pumping tends to produce precisely the situation where 
j3 = a 2 a 2 /2, and so it is predicted that MHD dynamical pumping will always cause configurations to tend 
towards being axially uniform (i.e., filamentary), hot, and dense, and with f3 = a 2 a 2 /2. 

The definition of a in Eq.(^) is closely related to that used for force- free fields. However, there is an 
important difference because while Eq.(p2|) corresponds to having p J z = aB z , Eq.(]22|) makes no statement 
about any relationship between and B^. The finite j3 equilibrium discussed in the previous two paragraphs 
is not force-free and involves the radial force balance J Z B$ = —dP/dr with = 0; this differs from the 
force-free radial equation J$B Z — J Z B^ = with /U J^ = aP^, p J z = «5 Z . It is worth noting that the 
determination of a made from vector magnetographs (e.g., Pevtsov et al., 1997) effectively use the definition 
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a = fi J z /B z which is equivalent to Eg. (p2|) ; these measurements do not provide information on either 
or Bj, and so do not provide any information on the value of \x J^lB^. Thus there is only one definition 
for a, but its application is different for force- free situations compared to finite-/? situations: for force- free 
situations a gives the ratio of current to flux for both toroidal and poloidal directions whereas for the finite 
(3 situation, a refers only to the the ratio of poloidal current to poloidal field. 

The prediction that (3 = a 2 a 2 /2 can be compared with the actual observed values of in TRACE 
flux loops. To calculate the predicted f3, we use the nominal measured flux loop radius a = 1.6 x 10 6 m from 
Aschwanden et al.(2000) and the nominal measured active region a = 2 x 1CP 8 m _1 from Fig. 4 of Pevtsov 
et al. (1997). These parameters give a nominal f3 pre dicted = ot 2 a 2 /2 = 5 x 10~ 4 . The observed value ^observed 
is calculated using a nominal measured density n = 10 15 m~ 3 and a nominal measured temperature 10 6 
K [Aschwanden et al. (2000)]. In addition a nominal axial magnetic field B z = 1.5 x 10~ 2 T is assumed 
based on the argument that because the flux tube is axially uniform, its axial field must also be axially 
uniform and so will have the same value as the nominal B z = 1.5 x 10~ 2 at the surface of an active region. 
These parameters give ^observed = l^onnT/B 2 = 4 x 10~ 4 which is very close to (3 pre dicted- If the model were 
wrong, one would expect no relationship between the predicted and observed /3's, i.e., one would expect a 
discrepancy of many orders of magnitude between the predicted and the observed /3's. 

This model also has implications regarding the brightening typically observed when the axis of a 
coronal loop starts to writhe and the loop develops a kink instability (sigmoid). Since kink instability 
occurs when ah ~ 1 and for a long thin flux tube a << h, this model predicts that (3 = a 2 a 2 /2 << a 2 h 2 /2 
will still be small even if a is increased to the point where ah ~ 1 and kink instability occurs. However, (3 
will increase as a increases and so this model predicts that the loop should brighten in proportion to the 
writhing of its axis (i.e., in proportion to a as ah approaches unity). 

Finally, we note that Feldman (2002) has recently used purely observational evidence to argue that electric 
currents with geometry similar to what is discussed here are the means by which the Sun and similar stars 
produce their coronal activity. 
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